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Prior modal stability analysis (Kojima et al. 1984) predicted that a rising or sedimenting 
droplet in a viscous fluid is stable in the presence of surface tension no matter how 
small, in contrast to experimental and numerical results. By performing a non-modal 
stability analysis, we demonstrate the potential for transient growth of the interfacial 
energy of a rising droplet in the limit of inertialess Stokes equations. The predicted 
critical capillary numbers for transient growth agree well with those for unstable shape 
evolution of droplets found in the direct numerical simulations of Koh & Leal (1989). 
Boundary integral simulations are used to delineate the critical amplitude of the most 
destabilizing perturbations. The critical amplitude is negatively correlated with the linear 
optimal energy growth, implying that the transient growth is responsible for reducing 
the necessary perturbation amplitude required to escape the basin of attraction of the 
spherical solution. 


1. Introduction 

The instability of capillary interfaces has long been an intriguing topic in fluid me¬ 
chanics. Perhaps one of the earliest investigated interfacial instability phenomena is the 
Rayleigh-Taylor instability, where a denser fluid located above a lighter one protrudes 
into the latter due to any arbitrary small perturbation of the initially flat interface. How¬ 
ever, this protrusion is not always observed when a droplet rises or sediments into another 
density-contrasted fluid. According to Hadamard (1911) and Rybzynski (1911), a spher¬ 
ical translating droplet is a solution of this problem in the Stokes regime, regardless of the 
presence or magnitude of the surface tension. What remains unknown, however, is the 
existence of other equilibrium shapes of the droplet and the influence of surface tension 
on the stability of the spherical solution. 

Experiments were conducted by Kojima et al. (1984) to examine this issue. Two pat¬ 
terns of shape instability were observed: depending on the viscosity ratio A, a protrusion 
or an indentation at the rear of droplet was seen to grow with time. Kojima et al. (1984) 
also performed a linear stability analysis assuming that the droplet underwent small de¬ 
formations. A linear operator depending on the viscosity ratio A and capillary number 
Ca (inversely scaling with the surface tension) was derived which governs the linearised 
droplet shape evolution. It was found that, irrespective of the value of Ca, i.e. even for 
arbitrary small surface tension, the eigenvalues of the operator had negative real part, 
pointing to a linearly stable shape. The authors recognized that this linear stability study 
contradicted their experiments showing instabilities with finite surface tension; Direct nu¬ 
merical simulations (DNS) (Koh & Leal 1989) also reported the unstable shape evolution 
of slightly disturbed droplets in the presence of sufficient surface tension (Ca < 10). Re¬ 
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cent numerical work has examined the effect of surfactants (Johnson & Borhan 2000) 
and viscoelasticity (Wu et al. 2012) on this scenario. 

The contradiction between the theory and experiments/DNS is somewhat reminis¬ 
cent of the case of the fingering instability of a film flowing down an inclined plane: 
the experimentally-measured (Huppert 1982; de Bruyn 1992) critical inclination angle 
triggering instability was found to be well below that obtained from the linear theory. 
Bertozzi & Brenner (1997) discovered that the traditional spectrum analysis failed to 
capture the short-time but significant energy amplification of the perturbations near the 
contact line. They pinpointed the missing mechanism by performing a so-called non- 
modal analysis, borrowed from the transient growth theory founded and developed in 
the early 1990s for hydrodynamic stability analysis (Trefethen et al. 1993; Reddy & 
Hcnningson 1993; Baggett et al. 1995), to identify and interpret the short-time energy 
amplification. 

The non-modal tools of stability theory have been used to explain the discrepancies 
between the theoretically computed critical Reynolds number and the experimentally- 
observed counterpart in a variety of wall-bounded shear flows (Schmid 2007). The tradi¬ 
tional eigenvalue analysis as also used in Kojima et al. (1984), i.e. the so-called modal 
approach, can sometimes fail to interpret real flow dynamics as the spectrum of the lin¬ 
ear operator only dictates the asymptotic fate of the perturbations without considering 
their short-term dynamics (Schmid & Henningson 2001). The non-modal analysis, in 
contrast, is able to capture the short-time perturbation characteristics and determine 
the most dangerous initial conditions leading to the optimal energy growth. In addition 
to its great success in the traditional hydrodynamic stability analysis, it has been also 
used to elucidate complex flow instability problems including capillary interfaces (Davis 
& Troian 2003), thermal-acoustic interactions (Balasubramanian & Sujith 2008; Juniper 
2011) and viscoelasticity (Jovanovic & Kumar 2010). 

In this paper, we perform a non-modal analysis to investigate the shape instability 
of a rising droplet in an ambient fluid, neglecting inertial effects. After introducing the 
linearised equations and operator in Sec. 2 and the non-modal approaches in Sec. 3, we 
demonstrate the existence of transient growth and predict the critical capillary numbers 
required for instability to become possible in Sec. 4. In Sec. 5, we conduct in-house DNS to 
compute the nonlinear shape evolutions of the droplets initiated with the linear optimal 
perturbations and identify the minimal amplitudes leading eventually to instability. We 
further analyse the relationship between the optimal growth and the critical amplitude 
of perturbation. We finally examine how the instability pattern is related to the viscosity 
ratio and propose a phenomenological explanation in Sec. 6 . 


2. Governing equations and linearisation 

We study the dynamics of a buoyant droplet rising in an ambient fluid in the Stokes 
regime. The droplet is assumed to be axisymmetric and the axis is along the z direction 
with gravity g = — ge,. The two Newtonian immiscible fluids, one carrying the droplet 
(fluid 2 ), and the other constituting the droplet (fluid 1 ) are characterized by different 
densities p 2 > pi, inducing (without loss of generality) an upward migration of the 
droplet. Likewise, their viscosities are p .2 and pi respectively, with a ratio A = pi/p 2 - The 
interface between the two fluids has a uniform and constant surface tension coefficient 7 . 
The undeformed state of the droplet is a sphere of radius a and terminal velocity t/t er = 
a g ^2~ Pl ^ 3(i+3A/2) (heal 2007). We use a and U ter as the reference length and velocity 
scales, and /^f/ter/a as the reference scale for p and <x, the modified pressure (removing 
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Figure 1: An axisymmetric droplet rising in the quiescent fluid, along the axial (z) di¬ 
rection. The fluid inside and outside is labelled as fluid 1 and fluid 2 respectively, so as 
their dynamic viscosities (pi, /r 2 ) and densities (p 1 , p 2 ). The polar coordinates R{8) are 
used to represent its shape, where 8 is measured from its rear stagnation point. L and B 
is the axis length along the revolution axis and orthogonal directions. 


the hydrostatic part) and the corresponding stress tensor respectively (Batchelor 2000). 
Hence, the governing equations for the non-dimensional velocity and pressure field inside 
the droplet (ui,pi) and that outside the droplet (u 2 ,p 2 ) are written as 

V • ui = 0, — Vpi + AV 2 Ui = 0, 

V • u 2 = 0,-Vp 2 + V 2 u 2 = 0, (2.1) 

where the velocity is zero at infinity and the boundary conditions on the interface are 


ui = u 2 , 

er 2 ■ n — ct\ • n = [V • n/Ca + 3z (1 + 3A/2) / (1 + A)] n. (2.2) 

Here, n is the unit normal vector pointing from the interface towards the carrier fluid 
and Ca = /r 2 [/ t er /7 is the capillary number indicating the ratio of the viscous effect with 
respect to the surface tension effect. 

Following Kojima et al. (1984), the interface of an axisymmetric droplet undergoing 
small deformation can be expressed in polar coordinates as 

E OO 

(2n + 1) f n P n (cos 8 ), (2.3) 

n=2 

where 8 is the polar angle measured from the rear of droplet, R ( 8) is the polar distance 
(see figure 1), S indicates the amplitude of the deformation, the P n are the nth-order 
Legendre polynomials and the f n are the corresponding coefficients. The first two terms 
Po and Pi are removed such that the volume of the droplet is conserved and its centroid 
stays at the origin (Kojima et al. 1984). To advance the interface, the kinematic condition 
dR(6 1 t)/dt = u (8,t) ■ n is applied. 

Following Kojima et al. (1984) and linearising the governing equations and truncating 
the series expansion, the evolution of the droplet can be obtained by solving a system of 
ordinary differential equations, 

df/d t = Af, (2.4) 

where the shape coefficient f = (/ 2 , fs ,,..., / m +i) T is a truncated vector and A is an 
m x m matrix depending on A and Ca. It should be noted that the shape of the droplet 
can be expressed by a unique series of coefficients f<5 and vice versa; for a certain f, the 
effective shape varies significantly with the amplitude and sign of 5. For the truncation 
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of f we use m = 1000 throughout our study: extensive tests using larger values of m 
confirm that our results are independent of this truncation level. 


3. Non-modal analysis: Theory 

As shown by the modal analysis of Kojima et al. (1984), the operator A has a sta¬ 
ble spectrum with all of its eigenvalues having negative real parts, irrespective of the 
magnitude of the surface tension, as long as the capillary number is finite. This model 
analysis predicts the long-term behaviour of the disturbance but in the short-term limit 
it is only valid if the linear operator A is normal, i.e. its eigenvectors are orthogonal. In 
the case of a non-normal operator, even though the amplitudes of all eigenmodes decay 
exponentially, their nonorthogonality can lead to a transient energy growth over a short 
time. We indeed found that A was non-normal. The optimal growth G max of the initial 
energy ( L 2 norm) over a chosen time interval [0,T] (Schmid 2007) is 


G max = max 

f(0) 


G (T) 


fcnik ' 

f(o)lk 


exp (TA)|| 2 , 


(3.1) 


where f (0) denotes the initial perturbation. G max ( T ) represents the maximum amplifica¬ 
tion of the initial energy at a target time (the so-called horizon) T where the optimization 
has been performed over all possible perturbations f (0). The optimal initial perturbation 
for horizon T will be denoted f™* (0). The quantity G max is the envelope of all individual 
gain profiles, indicating the presence of transient growth when G max ( T ) > 1 for some T. 

Compared with the L 2 norm in equ. 3.1, it is natural to introduce a physically-driven 
form of energy, designed for the physical problem at hand. In the present study, the vari¬ 
ation of surface area of the droplet AS is chosen as the target energy, since 7 A S indicates 
the interfacial energy throughout the evolution: AS is zero only for a spherical droplet 

and is positive otherwise. The surface area is S = 2n JJ R 2 sin 9\J 1 + [(\/R){dR/dO/f'dd. 
Assuming small deformation and thus a Taylor expansion yields 

-k^ 4 +k(f)> < 3 - 2) 


Plugging 2.3 into 3.2, the area variation AS = S — 47 t is found to be 

AS/ ( 2 ttS 2 ) = f T M AS f + o (<5 2 ) , (3.3) 

where Mas is the so-called weight matrix (Schmid 2001) of size m x m, with entries 

M A s (i, j) = 26^ (2i + 1) + ^ (2i + 1) (2j + 1) f P[ (cos0) P' (cost?) sin 3 ddd. (3.4) 

4 Jo 

The optimal growth of AS can now be defined as 

Vf T M A sf 


Gas x (T) = max 
v ' f(0) 


Gas (T) = 


VAS^rj 
y/ASl 0 ) 


= max 
f(0) 


Gas(T)= 


kf T (0)MAsf(0) 


(3.5) 


By Cholesky decomposition M A s = F r F, the above equation is formulated as 


G£s X (T) = max 
v 1 f(0) 


Gas (T) = 


|Ff (T) | 
11Ff (0) | 


(3.6) 
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In a similar way to how the asymptotic stability (t —> oo) is determined by the eigenval¬ 
ues of the evolution operator A, the maximum instantaneous growth rate of the perturba¬ 
tion energy at t = 0 + can be determined algebraically, expanding the matrix exponential 
exp(f A) r; I + tA at t = 0 + . The growth rate of the excess area AS is then 


1 dAS f T (0) [A t F t F + F t FA] f(0) 

A S~dT t=0+ “ f T (0)F T Ff(0) ' 

By introducing h = Ff(0), the maximum growth rate of AS is formulated as 


max 


1 dAS 


AS dt 


h T 


= max - 


FAF 


-l 


(FAF- 1 )" 




(3.8) 


which becomes the optimization of a Rayleigh quotient with respect to h. Because 
FAF” 1 + (FAF- 1 ) is a symmetric operator, the maximum is given by its largest eigen¬ 
value, 


max 


1 d^AS 
VAS dt 


'-’max 

- (faf- 1 + (FAF- 1 ) 1 ’) 

£—0+ 

2 V v ) 


(3.9) 


where s max [•] denotes the largest eigenvalue. This maximum instantaneous growth rate 
is commonly called the numerical abscissa (Trefethen & Embree 2005), which is closely 
linked to the numerical range Was (A, F) defined as the set of all Rayleigh quotients, 

Was (A, F) = {* : a = (FAF _1 p, p) / (p, p)}. (3.10) 


The numerical range is the convex hull of the spectrum for a normal operator (and is 
therefore always in the stable half plane z r < 0 for a stable operator) , but can extend 
significantly to even protrude into the unstable half-plane z r > 0 for stable non-normal 
operators. Its maximum protrusion is equal to the numerical abscissa and thus determines 
the maximum energy growth rate at t = 0 + . 


4. Non-modal analysis: results 

4.1. Transient growth and numerical range 

In figure 2, we show the optimal growth of the interfacial energy (T) for viscosity 
of ratios A = 0.5 and 5, varying the capillary number Ca. The threshold value of Ca to 
yield transient growth is between 4 and 5, in accordance with the rightmost boundary of 
the numerical range (see inset) depicted in the complex plane ( z r ,Zi ). The boundary is 
almost tangent to z r = 0 at Ca « 4.9 for A = 0.5 and Ca ~ 4.53 for A = 5 representing 
the critical capillary number Ca c above which the maximum energy growth rate at t = 0, 
maxf(o) y=^p| t= o+, is positive, guaranteeing transient growth. 

4.2. Linear growth and shape evolutions 

Non-modal analysis not only predicts the maximum energy growth over a particular time 
interval, but also provides the optimal perturbation, i.e. the initial shape coefficients 
f™ 4 (0) that ensure the optimal gain at horizon T. Figure 3 depicts the individual energy 

gains Gas for four optimal initial conditions fSS* (0) corresponding to T = 0.2, 1.05, 3.95 
and 5.45, with A = 0.5 and Ca = 6. Their gain profiles are tangent to G^cf (T) at t = T. 
The optimal perturbation targeting T = T max = 2.95 coincides with the optimal growth 
G^| x at its peak. 

Assuming small deformation amplitude and integrating equ. 2.4 in time, the linear 
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Figure 2: The optimal growth of the interfacial energy G^s x (T) versus the nondimen- 
sional time T, for viscosity ratio A = 0.5 (a) and A = 5 (b); for each case, four capillary 
number Cas are shown and for the highest Ca, the time T max corresponding to the peak 
energy growth is marked. The inset shows the boundary of the numerical range (z r , z%). 


6 < 0 T=T max =2.95 



S > 0 
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Figure 3: Linear growth Gas of the interfacial energy of the droplets with an optimal 
initial perturbation f™* (0) for the target times T = 0.2, 1.05, 2.95 and 5.45; the solid 
curve indicates the optimal growth G^| x (T) and it reaches its peak at T = T max = 2.95. 
The linear shape evolution of the perturbations are shown for negative and positive 5, 
on the left and right panel respectively. 


shape evolution is readily reconstructed for the droplets with the four optimal initial 
conditions, depicted in figure 3, at time t = 0 (dashed), 2.5,5,7.5,10 (light solid) and 
the target time t = T (solid); the evolution is shown for negative/positive <5 in (a)/(c). 
For both signs, the initial perturbation is mainly introduced near the tail (9 = 0) of the 
droplet where the interface is respectively flattened for S < 0 and stretched for S > 0 while 
the front part of the droplet remains spherical. In accordance with the modal analysis 
implying a linearly stable evolution, the perturbations eventually decay and the droplets 
finally recover a spherical shape. 


5. Nonlinear analysis 
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Figure 4: (a): Nonlinear energy growth Gas of the droplets with the optimal perturbation 
fpf, 4 (0) for the target time T = T max = 2.95; the solid curve indicates the linear energy 
growth. For positive 6 , Gas of droplets with <5 = 0.0496 and 0.0505 are shown, the 
former/latter being stable/unstable; for negative <5, the chosen value leading to stable 
and unstable evolution is <5 = —0.122 and S = —0.126 respectively, (b): The shape 
evolutions of the corresponding droplets. 


5.1. Nonlinear energy growth and shape evolution using DNS 

As the droplets deform more and more on increasing the initial perturbation ampli¬ 
tude, nonlinearities become significant and the droplet evolution cannot be adequately 
described by the linearised equations. We resort to DNS to address the non-linear dynam¬ 
ics using a three-dimensional axisymmetric boundary integral implementation, following 
the standard approach of Koh & Leal (1989). 

We focus on the droplets of A = 0.5 and Ca = 6 with the optimal perturbation fSF* , (0) 
achieving the peak of the optimal energy growth G^f* at T max . Two slightly different 
magnitudes of perturbation S = 0.0496, 0.0505 are chosen for the positive S and similarly 
S = -0.122,-0.126 for the negative case. Their energy growth Gas O = \J as(o) 
plotted in figure 4(a), together with the linear counterpart Gas if ) using equ. 3.3. The 
linear and non-linear energy growth share the same trend in the initial growing phase 
t < 3, but differ as the former is approximated by a truncated Taylor expansion. For the 
two values of <5 with the same sign but slightly different magnitude, the energy growth 
curves almost collapse before reaching their peaks at f ~ 4, but diverge afterwards; Gas 
decays for the smaller magnitudes S = —0.122 and <5 = 0.0496 indicating stable evolutions 
but maintains a sustained value around 1 for larger initial amplitudes S = —0.126 and 
S = 0.0505, implying the onset of instability. 

The shape evolutions of droplets are shown in figure 4(b). For <5 == —0.122 and —0.126, 
no significant difference is observed for 0 < t < 7.5: an inward cavity develops at the rear 
and sharpens; it is subsequently smoothed out and disappears for S = —0.122 while it 
keeps growing to form a long indentation for S = —0.126. These two values of <5 bound 
a threshold initial amplitude required to excite nonlinear instabilities. A similar trend 
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(b) (A, Ca) = (5, 6) 


Figure 5: The critical perturbation magnitude S c for: (a): (A, Ca) = (0.5,6) and (b): 
(A, Ca) = (5, 6). The upper and lower limits of S c (measured by the left scale) are plotted 
versus the target time T, with a curve fitted to show the trend. Accordingly, the linear 
energy growth Gas (measured by the right scale) is provided. 5^ and is the critical 
magnitude for an initially prolate and oblate respectively. 


is found for positive values of <5, while the instability arises through the formation of a 
dripping tail. 

It becomes natural to introduce <5 C , the critical magnitude of the perturbation above/below 
which the evolution of the drop is unstable/stable. Parametric computations are con¬ 
ducted to identify <5^ within a confidence interval (for instance 6+ £ [0.0496, 0.0505] and 
S~ £ [-0.122,-0.126] as in figure 4(a)). Searching in both directions, the critical am¬ 
plitude is then defined as S c = min(|£+|, |<$“|). When A = 0.5 and Ca = 6, |<5“| > |<5+|, 
implying that the instability tends to favour an initially stretched tail with respect to a 
flattened bottom; otherwise when A = 5 the situation reverses (|<5“| < | d+1), as discussed 
in next section. 


5.2. Critical amplitude of the perturbation 5 C 
Following the description of the previous paragraph, the critical deformation amplitude 
S c ( T ) can be determined for any targeting time T and associated optimal initial pertur¬ 
bation fp?* (0). The critical deformation amplitude S c ( T ) is plotted in figure 5 for Ca = 6 
and both A = 0.5 and A = 5, together with the optimal growth G^| x . The critical defor¬ 
mation amplitude S c is negatively correlated with G^cf and corresponds to a target time 
T slightly larger than T max where the peak transient growth is reached. This shows that 
the transient growth reduces the threshold non-linearity needed to trigger instabilities 
and consequently the critical magnitude of the initial perturbation. 

We also determined the critical amplitude Sc^° for an initially prolate (P) / oblate 
(O) ellipsoidal droplet to be unstable, as reported in figure 5. When the fluid inside the 
droplet is less viscous than the one outside, i.e. A < 1, an initially prolate droplet is more 
unstable, is less than half that of an oblate droplet S°; the trend reverses as A > 1. 
Such an observation is in agreement with the results of Koh & Leal (1989) using DNS 
(see fig. 11 of their paper). As expected, the minimum S c using the optimal perturbations 
is smaller than min((5^, <5(p) based on the limited family of ellipsoidal shapes. 

So far, we have analysed the critical amplitude S c of perturbations exhibiting transient 
energy growth. We would like to know how it varies as the transient growth decreases and 
even disappears as it is suppressed by high surface tension. In addition to Ca = 6, the 
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Figure 6: Akin to figure 5, adding 5 C of two smaller Cas for (a): A = 0.5 and (6): A = 5. 



Figure 7: The flow field co-moving with the droplet, using the optimal initial coefficient 
fpF j (0) , when (A, Ca) = (0.5,6), (a): <5 = —0.126 and (6): <5 = 0.0505. The red dot 
indicates the stagnation point of the flow. 


time-dependence of <5 C is shown in figure 6 for Ca = 4, 5. As expected, S c increases with 
decreasing Ca, by a factor of approximately 3, varying from the highest to the lowest Ca. 
With respect to T, S c varies non-monotonically for Ca = 5, 6 showing transient growth. 
In the absence of transient growth, like for Ca = 4, 5 C increases with T monotonically. 
Indeed, without transient growth, the energy decays monotonically and T max = 0, hence 
the minimum <5 C appears at T « 0. 


6. Conclusion and discussions 

In this paper, we have performed non-modal analysis and DNS to investigate the shape 
instabilities of an intertialess rising droplet which tends to recover the spherical shape, the 
attractor solution, due to surface tension. For sufficiently low surface tension, transient 
growth of the interfacial energy arises and leads to a bypass transition. This reduces the 
initial disturbance amplitude required to trigger instability, hence significantly decreasing 
the threshold magnitude of perturbation for the droplet to escape the basin of attraction. 
This magnitude is negatively correlated to the optimal growth of the interfacial energy. 
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We now compare our results with the work of Koh & Leal (1989) who employed 
DNS to identify the critical capillary number Ca c leading to shape instabilities of an 
initially prolate or oblate ellipsoidal sedimenting droplet; the magnitude of perturbation 
is A = (see figure 1). For their lowest magnitude | A| = 1/21 considered, Ca c € (4, 5) 
for A = 0.1,0.5 and 5, indeed close to our prediction: Ca c ~ 5.42, 4.9 and 4.53 respectively 
for the same A. Additionally, Koh & Leal (1989) observed that for a viscosity ratio 
A < 1/A > 1, the first unstable pattern appears as a protrusion/indentation developing 
near the tail of the droplet that is initially a prolate/oblate. The trend holds in our case 
even though we search over all possibilities for the most ’dangerous’ initial perturbation 
instead of using an initially ellipsoidal shape. This is also reflected from the initial shapes: 
as A < 1/A > 1, the optimal shape shares a common feature with an oblate/prolate 
ellipsoid, namely its rear interface is compressed/stretched. 

To explain the dependence of the instability patterns on the viscosity ratio A, let us 
focus on the velocity field near the tail of the droplet (see figure 7), where the flow 
resembles a uniaxial extensional flow, drawing the tip into the drop on the top side and 
pulling it outwards on the other side. We suggest that this imbalance induces the onset 
of the shape instability. The internal (respectively external) viscous force on the tip is 
Hidu^v/dz (respectively /dz). When /zi < /z 2 , i.e. A < 1, the external viscous 

effect overcomes the internal one, hence the perturbation tends to be stretched outward 
to form a protrusion; otherwise, when A > 1, it is prone to be sucked inwards to form an 
indentation. 

Developed originally for hydrodynamic stability analysis, non-modal tools have here 
demonstrated the predictive capacity for the inertialess shape instabilities of capillary 
interfaces. This work might stimulate the application of non-modal analysis for complex 
multiphase flow instabilities even at low Reynolds number. 


Acknowledgements 

L.Z. thanks Francesco Viola for the helpful discussions. We thank the anonymous 
referee for pointing out an incorrect coefficient in our previous derivation. Computer 
time from SCITAS at EPFL is acknowledged, and the European Research Council is 
acknowledged for funding the work through a starting grant (ERC SimCoMiCs 280117). 


REFERENCES 

Baggett, J. S., Driscoll, T. A. & Trefethen, L. N. 1995 A mostly linear model of tran¬ 
sition to turbulence. Phys. Fluids 7 (4), 833-838. 

Balasubramanian, K. & Sujith, R. I. 2008 Thermoacoustic instability in a rijke tube: Non¬ 
normality and nonlinearity. Phys. Fluids 20 (4), 044103. 

Batchelor, George Keith 2000 An introduction to fluid dynamics. Cambridge university 
press. 

Bertozzi, A. L. & Brenner, M. P. 1997 Linear stability and transient growth in driven 
contact lines. Phys. Fluids 9 (3), 530-539. 

de Bruyn, J. R. 1992 Growth of fingers at a driven three-phase contact line. Phys. Rev. A 
46 (8), R4500. 

Davis, J. M. & Troian, S. M. 2003 On a generalized approach to the linear stability of spatially 
nonuniform thin film flows. Phys. Fluids 15 (5), 1344-1347. 

Hadamard, J. 1911 Mouvement permanent lent d’une sphere liquide et visqueuse dans un 
liquide visqueux. C. R. Acad. Sci. Paris 152, 1735-1738. 

Huppert, H. E. 1982 Flow and instability of a viscous current down a slope. Nature 300 (5891), 
427-429. 



Inertialess non-modal energy growth of a rising droplet 


11 


Johnson, Robert A & Borhan, Ali 2000 Stability of the shape of a surfactant-laden drop 
translating at low reynolds number. Physics of Fluids (1994-present) 12 (4), 773-784. 

Jovanovic, M. R. & Kumar, S. 2010 Transient growth without inertia. Phys. Fluids 22 (2), 
023101. 

Juniper, M. P. 2011 Triggering in the horizontal rijke tube: non-normality, transient growth 
and bypass transition. J. Fluid Mech. 667, 272-308. 

Koh, C. J. & Leal, L. G. 1989 The stability of drop shapes for translation at zero reynolds 
number through a quiescent fluid. Phys. Fluids 1 (8), 1309-1313. 

Kojima, M., Hinch, E.J. & Acrivos, A. 1984 The formation and expansion of a toroidal drop 
moving in a viscous fluid. Phys. Fluids 27 (1), 19-32. 

Leal, L. G. 2007 Advanced transport phenomena: fluid mechanics and convective transport 
processes. Cambridge University Press. 

Reddy, S. C. & Henningson, D. S. 1993 Energy growth in viscous channel flows. J. Fluid 
Mech. 252, 209-238. 

Rybzynski, W. 1911 Uber die fortschreitende Bewegung einer fliissingen Kugel in einem zahen 
Medium. Bull. Int. Acad. Sci. Cracovie 1A, 40-46. 

Schmid, P. J. 2001 Tools for nonmodal stability analysis. Notes from a Ladhyx Tutorial p. 24. 

Schmid, P. J. 2007 Nonmodal stability theory. Annu. Rev. Fluid Mech. 39, 129-162. 

Schmid, P. J. & Henningson, D. S. 2001 Stability and transition in shear flows, Applied 
Mathematical Sciences, vol. 142. New York: Springer-Verlag. 

Trefethen, L. N. & Embree, M. 2005 Spectra and pseudospectra: the behavior of nonnormal 
matrices and operators. Princeton University Press. 

Trefethen, L. N., Trefethen, A. E., Reddy, S. C. & Driscoll, T. A. 1993 Hydrodynamic 
stability without eigenvalues. Science 261 (5121), 578-584. 

Wu, H, H A j-Hariri, H & Borhan, A 2012 Stability of the shape of a translating viscoelastic 
drop at low reynolds number. Physics of Fluids (1994-present) 24 (11), 113101. 


